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Abstract 

What range of momentum components in the deuteron wave function are 

available ed elastic scattering data sensitive to ? This question is addressed 

within the context of a model calculation of the deuteron form factors, based 

on realistic interactions and currents. It is shown that the data on the A{q), 

B{q), and T2o(g) observables at g < 6 fm~^ essentially probe momentum 

components up to ~ Am.,^. 
21.30.+y, 21.45.+V, 27.10.+h 
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I. INTRODUCTION, RESULTS, AND CONCLUSIONS 



The present work addresses the following issue: what range of momentum components 
in the deuteron wave function are probed by presently available e-d elastic scattering data ? 
The question is answered within the context of a model [|l|, based on realistic interactions and 
currents, that predicts quite well the observed deuteron A{q) and B{q) structure functions 
and T2o{q) tensor polarization up to momentum transfers g ~ 6 fm^^. To this end, the 
deuteron S- and D-wave components obtained in the full theory and denoted as ul{p), 
L=0, 2, are truncated, in momentum space, as 

UL{p;n) = - , , ,1/2 Ul{p) , (1-1) 

[1 + exp{p - pn)/a\ 

where the cutoff momentum pn=nmT^ (m^r is the pion mass and n is an integer) and 
a=0.1m^. The constants c„ are fixed through the normalization condition, Eq. ( p.6| ) below. 
The wave functions UL{p',n) corresponding to n=l, 2 and 4 along with the reference ul{p) 
are shown in Fig. |1|. The A{q), B{q), and T2o{q) observables are calculated using these 
truncated wave functions by including one-body only and both one- and two-body currents. 
The results are displayed in Figs. 

Aspects of the present theory [|I| pertaining to the input potential, the boost corrections 
in the wave functions, the one and two-body charge and current operators, are succinctly 
summarized in Sees. |I^^V|. Here it suffices to note that the results shown in Figs. are 
obtained within a scheme in which these various facets are consistently treated to order 
(f/c)^. Thus the present approach improves and extends that adopted in Refs. and 0, 
in which relativistic corrections of order (f/c)^ were selectively retained. For example, 
the (f/c)^ contributions associated with the vr-exchange charge operator were included, 
while those originating from the one-pion-exchange potential and from boosting the wave 
function were ignored. The consistent (f/c)^ scheme employed here, however, does not alter 
quantitatively the predictions obtained in Refs. except for the B{q) structure function 
as discussed in Sec. |V|. A detailed analysis of the differences between the two approaches is 
beyond the scope of this work; it will be presented in Ref. 



We assume that the three-momentum transfer q in e-d elastic scattering is in the z- 
direction. In the Breit frame the deuteron has initial and final momenta — g/2 and +q/2 
in this direction, respectively. In the presence of only single-nucleon currents, the deuteron 
wave function has to have components with relative momenta p=\ pi — P2 | /2 larger than 
g/4 in order to produce elastic scattering at momentum transfer q. In a limiting case, for 
example, the nucleons have momenta of —q/2 and in the initial state and +q/2 and in the 
final, all along the z-axis. Thus elastic scattering via one-body currents probes momentum 
distributions ul{p) a.t p > g/4. This argument, however, does not establish the maximum 
value of p which the observed form factors with g < 6 fm~^ are sensitive to. In addition, 
there is no kinematical limit on initial and final relative momenta for scattering via pair 
currents. Therefore, we use a realistic model of the deuteron to estimate the maximum 
value of p probed by the available data. 

There is growing interest in chiral effective field theories of nuclear forces and cur- 
rents @,|^. In such theories the effective Lagrangian is obtained after integrating out states 
with momenta greater than a specified cutoff. We hope that the present results provide an 
estimate of the deuteron structure obtained in theories with cutoff momentum pn=nmT^. 

The deuteron wave function is dominated by the one-pion exchange potential f7r(k), 
where k=p' — p is the momentum transferred by the exchanged pion. The dependence of 
f7r(k) (see Eq. ( |2.2|) below) on the magnitude k is primarily given by the ttNN form factor 
fn{k) for k > itLt^. In this limit the leading non-relativistic term of f7r(k) becomes: 



expected to be sensitive to f-^^k) up to A; ^ p. 

The results shown in Figs. indicate that the available data on deuteron form factors 
at g < 6 fm^^ confirm the present wave functions up to p ~ 4m7r. After including pair- 
current contributions, the reference wave function and that truncated at p=4m^ do not 




f^^2(fc)^i.ka2-kTi-r2, 




where k is a unit vector. Thus the deuteron wave function at relative momentum p is 
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give significantly different form factors in tliis g-range. Tlie reference and truncated wave 
functions, respectively, over- and under-predict tlie observed A{q) at g > 6 fm~^; botli 
under-predict B{q < 7 fm~^), and fail to reproduce tlie data point for T20 at q 6.8 fm^^. 
The available data on A{q) at larger values of q can be used to test the deuteron wave 
function aX p > 4:171^^; however, improved theoretical understanding and more accurate data 
on B{q > 6 fm~^), and new data on T2o{q > 6 fm~^) are needed. 



II. DEUTERON WAVE FUNCTION 

The deuteron rest-frame wave function is obtained by solving the momentum-space 
Schrodinger equation with the relativistic Hamiltonian 



H^" = 2^p^ + m? + , (2.1) 

where consists of a short-range part parameterized as in the Argonne fig potential 0, 
and of a relativistic one-pion-exchange potential (OPEP) given by 



<(p;.Pi)= 



CTi ■ kcr2 ■ k 

T1-T2. (2.2) 



+ k"^ E' E 



E' + m E + m 

Here m is the nucleon mass, fnNN is the pion-nucleon coupling constant (/^^^/47r=0.075), 
Pi and p[ are the initial and final momenta of particle 1 in the center-of-mass frame, k=p'i — 



Pi is the momentum transfer, E=ypi + m^, and E'=yp[ + m^. The monopole form factor: 

/\2 _ ™2 

with A^=1.2 GeV/c is used in this work. 

The /i-dependent term characterizes possible off-energy-shell extensions of OPEP, and 
leads to strong non- localities in configuration space. In particular, the value fi=-l (/^=1) 
is predicted by pseudoscalar (pseudovector) coupling of pions to nucleons, while fi=0 cor- 
responds to the so-called "minimal non-locality" choice f^. It has been known for over two 



decades 0, and recently re-emphasized by Forest that these various off-shell extensions 
of OPEP are related to each other by a unitary transformation, in the sense that 



1/i 



H 



U 



(2.4) 



if terms of 27r-range (and shorter-range) are neglected. The hermitian operator U is given 
explicitly in Ref. 0. This unitary equivalence implies that predictions for electromagnetic 
observables, such as the form factors under consideration here, are independent of the par- 
ticular off-shell extension adopted for OPEP, provided that the electromagnetic current op- 
erator, specifically its two-body components associated with pion exchange, are constructed 
consistently with this off-shell parameter. This point will be further elaborated below. At 
this stage, however, it is important to recall that Forest has constructed relativistic Hamil- 
tonians H^^ with /i=±l, 0, each designed to be phase-equivalent to the non-relativistic H, 
based on the Argonne fig potential. 

Given this premise, the /i=0 prescription is adopted for OPEP from now on, and the 
superscript is dropped from H^^^, for ease of presentation. The resulting momentum-space 
wave function in the rest frame of the deuteron is denoted with V'Af (p; 0) (here, the argument 
indicates the rest frame in which the deuteron has velocity v=0), and is written as 



ijMip; 0) = %(p)3^on(p) + «2(p)3^2n(P 



(2.5) 



where p is the relative momentum, M is the angular-momentum projection along the z- 
axis, tjq is the pair isospin T=0 state, and y^sjip) standard spin-angle functions. The 
normalization is given by: 



dpp^ 



Hip) + ulip)] = 1 



(2.6) 



/o (27r)3' 

The internal wave function in a frame moving with velocity v with respect to the rest 
frame can be written to order (f/c)^ as 



^m(p; v) 



1-^ 



1 



1 - ^(v ■ p)(v ■ Vp) - ■ (<Ti - 0-2) X p 



V^m(p;0) 



1 - —V ■ (CTI - 0-2) X p 

4m 



V'A/(Pii/7,p±;0) 



(2.7) 



where 7 = — f 2, and py and denote the components of the momentum p parallel 

and perpendicular to v, respectively. Only the kinematical boost corrections, associated 
with the Lorentz contraction (term oc v ■ pv ■ Vp) and Thomas precession of the spins 
(term oc v ■ (cri — 02) x p) are retained in Eq. (|2.7] ). There are in principle additional, 
interaction-dependent boost corrections. Those originating from the dominant one-pion- 
exchange component of the interaction have been constructed explicitly in Ref. p| to order 
[vjcf'. At this order, however, their contribution to the form factors in Eqs. (f4.5|) - (|4.7|) 
vanishes, since it involves matrix elements of an odd operator under spin exchange between 
S=\ states. Indeed, this same selection rule also holds for the Thomas precession term. 
Therefore the interaction-dependent boost corrections do not contribute to elastic scattering 
up to order (t;/c)^, and are neglected in the following. 

It should also be noted that terms of order higher than (f/c)^ due to Lorentz contraction 
have in fact been included in the second line of Eq. ( |2.7| ). The factor 1/^^ ensures that the 
wave function in the moving frame is normalized to one, ignoring corrections of order [vjc]'^ 
from the Thomas precession term. 

It is interesting to study the relation between the Breit-frame matrix element p(q; -B) of 
the point-nucleon density operator, and the rest frame p(q;0). One finds, again ignoring 
Thomas precession contributions: 

1 /■ dp . { pii + q/4 \ f pii - q/4 \ 

= PM(q/7B;0) , (2.8) 

after rescaling of the integration variables, (p||/7b,P±) — ^ (PihPx), in the last integral. 
Here 7^ is the Lorentz factor corresponding to = (g/2)z/-y/g2/4 + in the Breit frame. 
This result is consistent with the naive expectation that the density in configuration space 
is "squeezed" in the direction of motion by the Lorentz factor 7^ or, equivalently, that 
its Fourier transform is "pushed out" by 75. For g=6 fm"^, vb — 0.3 and 75 ^ 1.05, 
corresponding to a 5 % Lorentz contraction. 
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III. NUCLEAR ELECTROMAGNETIC CURRENT 



The electromagnetic charge and current operators include one- and two-body terms. 
Only the isoscalar parts need to be considered, since the deuteron has T=0. The one-body 
term is taken as 



ir-body(q) = (Pl-body(q),jl-body(q)) = ^ j 



i=l,2 2 



Ui , (3.1) 



where g'^=(0,gz) in the Breit frame, Ui and are the initial and final spinors of nucleon 

i, respectively, with Ui = m|7°, and and F2 denote the isoscalar combinations of the 

nucleon's Dirac and Pauli form factors, normalized as Ff{0)=l and -F2^(0)=— 0.12 (in units 

of n.m.). The Hohler parameterization of Fi and F2 is used in this work. The spinor u, 

or rather its transpose, is given by 

E + m\^/V cr-p 



where p and E=y/p^ + rn? are the nucleon's momentum and energy, and Xs is its (two- 
component) spin state. Note that u^u=l. In earlier published work on the form factors 
of the deuteron |0J^ and A=3-6 nuclei, most recently [|l^,|ri|], boost contributions were 
neglected, and only terms up to order (f /c)^ were included in the non-relativistic expansion 
of 

jT-body5 namely the well known Darwin- Foldy and spin-orbit corrections to p(q). In the 
calculations reported here, however, the full Lorentz structure of j^-body retained. 

The two-body terms included in j2-body those associated with vr- and p-meson ex- 
changes and the p7r7 transition mechanism. The vr and p spatial current operators to leading 
order in f/c are isovector, and therefore their contributions vanish in e-d elastic scattering. 
The TT-exchange charge operator is obtained, consistently with the off-shell parameter ix 
adopted for OPEP, as [| 

= ^^^^fT^"^ • q .2 ■ k^r, . r2 + 1 ^ 2 , (3.3) 

where kj=p- — pj, /^(A;) is the vrA^A^ form factor defined in Sec. 0, and Ff{q) denotes the 
Hohler parameterization for the nucleon's isoscalar Dirac form factor. It is easy to show, 
to order (f/c)^, that 



-ifiU 



(Pl-body + ^'^^ - Pl-body + + iP [Pl-body , ^] ^ Pl-body + P^ , (3.4) 



where exp(— iyuf/) is the unitary transformation of Sec. ||. Again, we emphasize that the 
choice /i=0 is made for (as for OPEP in Sec. 0) in the present work. 

At this point it is important to recall that the vr-exchange charge operator corresponding 
to p=-l was included in all our previous studies of light nuclei form factors (for a review, see 
Ref. []T2[), based on non-relativistic Hamiltonians, in which the OPEP part of v was taken 
to be its leading local form, i.e. E'=E and E/m ^ 1 in Eq. While these calculations 

are not strictly consistent, since relativistic corrections are selectively retained in p(q) but 
ignored in the potentials and wave functions, they nevertheless give fairly accurate results 
since the corrections to the wave function are rather small 

Vector-meson (p and uj) two-body charge operators have been found to give contributions. 



in calculations of light nuclei form factors [|T^, that are typically an order of magnitude 
smaller than those associated with the p^r operator. In the present study, only the p-exchange 
charge operator is considered: 

^'^"^^ " 8^;^^ ' ^'^ \{+w? ^'^^ X q) • (^^2 X ka) Ti ■ T2 + 1 ^ 2 , (3.5) 

where mp, QpNN and i^pNN are the p-meson mass, vector and tensor coupling constants, 
respectively. Here the values 5'p^^/(47r)=0.84, and KpMN=^-^ are used from the Bonn 2000 
potential while the pNN monopole form factor fp{k) has Ap=1.2 GeV/c. 
Finally, the pirj current is obtained from the associated Feynman amplitude as 

,M (a)-\n 3pNN UnN ^ / n fpjkl) fAh) p^arj 



X 



— / / I • f^pNN J a \ 



(^^2 7/3^2 75^i2) Ti ■ T2 + 1 ^ 2 , (3.6) 



where gpT,.y and GpT,^{q) are the pn'y coupling constant and form factor, respectively, and 
eoi23=l. The value (7p^^=0.56 is obtained from the measured width of the decay p ^ n'j [jl4|, 
while the form factor is modeled, using vector-meson dominance, as 

GpTT-riq) = — — T7 — 2 ' ^'^■'^^ 
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where rriuj is the w-meson mass. Note that in Eq. ( p.6|) the term proportional to fc"A;f in 
the p-meson propagator has been dropped, since it vanishes if the nucleons are assumed to 
be on mass shell. Furthermore, in both meson propagators retardation effects have been 
neglected. 

The leading terms in a non-relativistic expansion of Eq. ( |3.6| ) reduce to the familiar 
expressions for the pn'-f charge and current operators, as given, for example, in Ref. 0. 
However, it is known that these lowest-order expansions, particularly that for jp^^^, are in- 
accurate. This fact was first demonstrated by Hummel and Tjon [^] in the context of 
relativistic boson-exchange-model calculations of the deuteron form factors, based on the 
Blankenbecler-Sugar reduction of the Bethe-Salpeter equation. It was later confirmed, in a 
calculation of the deuteron B{q) structure function p|, that next-to-leading-order terms in 
the expansion for jp^^, proportional (1 -|- KpArAr)/m^, very substantially reduce the contribu- 
tion of the leading term. This issue will be returned to later in Sec. 0. Here, again we stress 
that the full Lorentz structure of jj^^ is retained in the present study. 

IV. DEUTERON FORM FACTORS 

The deuteron structure functions A and B, and tensor polarization T20 are expressed in 
terms of the charge {Gq and G2) and magnetic (Gi) form factors as 

A{q)=Gl{q) + lvGl{q) + lv'Gl{q) , (4.1) 

B{q) = ^7l{l + v)Gl{q) , (4.2) 
^ / N ^x(x + 2) +y/2 , , 

= -^TWT^ • <^-^' 

where r,=(5/2m,,)^ x=(2/3)nG2(g)/G„{g), t,=(2/3)r, [1/2 + (1 + r,) tg20/2] |Gi(9)/G„M]", 
and 6 is the electron scattering angle. The expressions above are in the Breit frame, and 
therefore q denotes the magnitude of the three-momentum transfer. The form factors are 
normalized as 

Go(0) = 1 , Gi(0) = (m,/m)/i, , ^2(0) = m^Q, , (4.4) 



where fid and Qd are the deuteron magnetic moment (in units of n.m.) and quadrupole 
moment, respectively, and are related to Breit-frame matrix elements of the nuclear electro- 
magnetic charge, p(q), and current, j(q), operators via: 

Go{q) = 1 {i^M^^B I p(q) I ipM] -vb) , (4.5) 

"J A/=±1,0 

Gi{q) = — ^(^Af=i; vb I jA=i(q) I i'M=o; -^b) , (4.6) 

G2{q) = — [(V'M=o; vb I p(q) | V^m=o; -vb) - {^m=i; vb | p(q) | i'M=i; -vb)] • (4.7) 
2r] 

Here jx=i denotes the standard +1 spherical component of the current operator. 

The calculations are carried out in momentum space. The wave functions iPm{p] ^^b) are 
given in the second line of Eq. ( |2.7| ), and the charge and current operators are those described 
in Sec. The one-body current matrix elements involve the evaluation of integrals of the 
type, in a schematic notation, 

/ (^^^'^'*^^ ^ ^^)^'i^-body(p'i' Pi) ^m(p - q/4; -Vb) , (4.8) 

with pi=p + q/2 and pi=p — q/2, which are performed by standard Gaussian integrations. 
The two-body current matrix elements, instead, require integrations of the type 

/ (2^)3 (2^)3 ^^^''-P'' ^'2-body(ki, ka) iIjm{p; -vb) , (4.9) 

with ki=q/2 + p' — p and k2=q/2 — p' + p. These integrations are efficiently done by Monte 
Carlo techniques by sampling configurations (p, p') according to the Metropolis algorithm 
with a probability density W{p,p')=w{p)w{p'), where w{p)=p'^ \ 4'm=o{p',0) p. The com- 
puter programs have been successfully tested by comparing, in a model calculation which 
ignored boost corrections and kept only the leading terms in the expansions for jf_body 
i2-body5 present results with those obtained with an earlier, configuration-space version 
of the code P]. 
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V. FURTHER RESULTS 



In this section we briefly discuss the contribution of the prey current to the B{q) structure 
function. In Fig. ^ the results calculated with the current given in Eq. ( |3.6D (curve labelled 
p7r7-R) are compared with those obtained by using the leading term in its non-relativistic 
expansion (curve labelled p7r7-NR), 

+ 1^2. (5.1) 

Figure ^ demonstrates the inadequacy of the approximation ( |5.1| ), a fact which, as mentioned 
already in Sec. [TTTI , has been known for some time [jl5| . Indeed a more careful analysis shows 
that the contributions of next-to-leading order terms are not obviously negligible, since they 
are proportional to the large pNN tensor coupling constant, Kp7VAr=6.1 in Ref. |l^. We 
sketch the derivation again here for the sake of completeness. 

The vector structure eiy^r^" occurring in Eq. (|3.6|) , with F'^ = (F°,r) defined as 

r = Ii;(7^ + i^a'^"fci„)«i, (5.2) 

can be written in the Breit frame as 

err.arT" k^ = F° (q X ki) - A:? (q X r) , (5.3) 

where the index i = 1, 2, 3 and q=ki + k2. Up to order (f/c)^ included, one finds, in analogy 
to the non-relativistic expansion of the nucleon electromagnetic current (with point-nucleon 
couplings, i.e. F^{q) 1 and F^{q) Kp^n) 

kl i 



V' -1-{1 + 2kpnn) 



o-i ■ IPi X Pi. 



(5.4) 



^ ~ 2m ^"^'^ ^^'^ 2m "^"^ '^'"^^'^ cri x ki . (5.5) 
Retaining only the leading term in F'^, namely F*^ ~ (1,0), leads to the expression in 
Eq. (|5.1|). Some of the corrections proportional to npNN were explicitly calculated in Ref. 0, 
and were found to decrease substantially those from the leading term. 
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We note, in passing, that the p7r7 charge operator, Pp7r7(q), is proportional to 

Q,,,r^fc^g" = r-ki X q 

^ ^ (1 + i^pNn) cTi ■ ki X ks , (5.6) 

where the small non-local term proportional to p'^ + pi has been neglected. The standard 
form of the Pp7r7(q) commonly used in studies of light nuclei form factors [10,11| easily follows. 



In this case too, however, higher order corrections included in j"^^, Eq. (|3.6|) , reduce the 
contribution of the leading term, although they do not change its sign. A more detailed 
discussion of this issue will be given in Ref. 

We conclude this section with a couple of remarks. The first is that the destructive 
interference between the one-body and p7r7 currents is also obtained in recent calculations 
of the B{q) structure function, carried out in the covariant framework based on the spectator 
equation ||T6|,p^. 



The second remark is that in all earlier studies of light nuclei form factors [12] additional 
two-body currents, originating from the momentum-dependent terms of the two-nucleon 
potential, were included. Work is in progress ||T|] on an improved treatment of these currents, 
in an approach similar to that proposed in Ref. 
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FIGURES 




P{fm"') 

FIG. 1. The truncated wave functions ULip]n), with n=l, 2 and 4, are compared with the 
reference ul{p). 




q(fm-') 

FIG. 2. The A{q) structure functions obtained with the truncated wave functions UL{p',n) for 
71=1,2, and 4 are compared with the reference A(q) and data. The results obtained with one-body 
only and both two- and two-body operators, labelled respectively lA and TOT, are shown. 



14 



10" 



-3 



10 

10" 

10-^ 
§ 10-^ 

00 

10" 

10-^ 



10" 



10 



-10 



o Mainz-81 
□ Saclay-85 

Bonn-85 
A SLAC-90 j 

lA 

TOT 




01 234.5678 



q(fm-') 



FIG. 3. Same as in Fig. but for the B(q) structure function. 
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FIG. 4. Same as in Fig. |2|, but for the T2o(g) tensor polarization. 
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FIG. 5. The B{q) structure functions obtained with the relativistic and non-relativistic forms 
the piT'y current of Eqs. ( |3.6| ) and (|5.1|), labelled respectively pirj-R and p7r7-NR. 
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